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Abstract 

We show, by explicit computation, that bare lattice perturbation the- 
ory in the two-dimensional 0{n) nonlinear a models with superinstan- 
ton boundary conditions is divergent in the limit of an infinite number 
of points |A|. This is the analogue of David's statement that renormal- 
ized perturbation theory of these models is infrared divergent in the 
limit where the physical size of the box tends to infinity. We also give 
arguments which support the validity of the bare perturbative expan- 
sion of short-distance quantities obtained by taking the limit |A| — > oo 



*Work supported in part by Schweizerischer Nationalfonds 

^On leave from the Institute of Theoretical Physics, Eotvos University, Budapest, 
Hungary 



term by term in the theory with more conventional boundary condi- 
tions such as Dirichlet, periodic, and free. 



2 



1 Introduction 



Quantum chromodynamics (QCD) is the presently favored candidate the- 
ory for strong interactions. However, to estabhsh its esteemed status the 
theory must be able to reproduce the low-lying spectrum and low-energy 
S-matrix elements to a respectable precision. To accomplish this we require 
a nonperturbative definition of the theory - the most promising of which 
is lattice regularization. On the other hand, many aspects of high energy 
phenomena involving hadrons (e.g., jets, deep inelastic scattering, etc.) are 
described successfully using renormalized perturbation theory (PT); the ra- 
tionale being the expected property of asymptotic freedom, i.e., that the 
amplitudes can be expressed as a series in a coupling which depends on the 
energy of the process and goes to zero as the energy goes to infinity. The 
conventional wisdom is that "in principle" one could justify this use of PT by 
showing that the coefficients of the conventional perturbative series provide 
the true asymptotic expansion of the full theory defined nonperturbatively 
via the lattice regularization. Confirming this proclaimed status of PT is 
also essential in aiming at a criterion where to truncate it. Without such 
a criterion PT produces sequences of predictions, one for each order where 
one decides to truncate it, and a "priori none of them might be close to the 
true answer. It must be emphasized, however, that the usual arguments for 
the above scenario are far off a proof and Patrascioiu and Seller repeatedly 
emphasized that also alternative scenarios can be imagined [||, 

In full QCD a rigorous understanding of the status of PT is likely not 
to become amenable in the near future, so that it is sensible to address the 
question in a simpler context, e.g., the 0(n) nonlinear a models in two di- 
mensions. These models are also perturbatively asymptotically free and are 
thought to describe a multiplet of stable massive particles, the mass scale 
being nonperturbative in the coupling constant. In the lattice formulation 
one usually starts with a finite number of lattice points |A| and the standard 
discretization (|2.l| ). The conventional picture is that the critical point at 
which the continuum limit for n > 2 should be taken is /?c = oo. Patrascioiu 
and Seller on the other hand, conjecture that there is a critical /?c < oo 
(for all n > 1) beyond which the theory is in a massless phase. It might then 
happen that although the continuum limit obtained by approaching (i^ from 
below has a mass gap and describes the expected low-energy properties, the 
theory thus obtained is not asymptotically free at high energies. In par- 
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ticular, conventional perturbation theory in the infinite volume limit would 
then be wrong for 2-d models with a non-Abelian global symmetry. Their 
arguments are based on results from percolation theory |Q] and with analogy 
to the case of the Abelian model n = 2 (for which it is rigorously proven 
that there is a phase transition at finite /? 1^). Numerical investigations can 
give useful hints but obviously cannot establish which scenario is correct. 

Concerning the status of perturbation theory it is known that for a fixed 
number of lattice points |A| ~ L? and for a large class of boundary conditions 
(BC's) observables have a well-defined perturbative expansion in the bare 
coupling I3~^ . We shall use the term "short-distance observable" to denote 
0(re)-invariant correlation functions of some local field, where (in the bulk 
of the paper) all arguments are fixed lattice sites, a distance 0{L) away from 
the boundary. The problems now arise with the limit L — > oo. Let us first 
consider the following question: 

Ql: Is the asymptotic expansion of short-distance observables for fixed L 
and /3 — > cxD uniform in L? Equivalently, can for such quantities the 
limits /3 — > oo and L — > oo be exchanged in bare lattice PT? 

A major finding of 0| is that the answer to Ql actually depends on the 
choice of boundary conditions. The reasoning can be illustrated with the 
computation of the energy density £■ at the center of the lattice. Although 
£■ does not have a continuum limit relevant to the physics, by application of 
the Mermin- Wagner theorem |p it should be independent of the boundary 
conditions in the limit of an infinite number of lattice points. Patrascioiu and 
Seller invented boundary conditions (which they coined superinstanton (SI) 
boundary conditions; c.f. sect. 2) and computed the one-loop coefficient of £. 
They showed that it has a finite limit as L ^ oo but the result was different 
from that with more conventional BC's such as Dirichlet. Furthermore they 
claim a similar result holds for the renormalization group /3-functions. The 
conclusion is that for at least one of the BC's involved the answer to Ql is 
negative. The authors interpret their result as casting doubt on the validity 
of standard PT for all BC's. They claim that the correct perturbative 
treatment in the limit of infinite volume must include expansions around 
the so-called superinstanton gas (configurations whose energy tends to zero 
as L — > oo; see sect. 3). 
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We believe that such conclusions are too strong and we set out to clarify 
some of the points raised. In the next section we show that the SI BC are 
"sick" in the sense that the limits in Ql cannot be interchanged. We do 
this by showing that the perturbative coefficients of E with SI BC diverge at 
two- loop order in the limit L — > c«.[| We have not considered the particular 
/3-functions defined in [Q] but we believe a similar result will be found. This 
infrared divergence of the PT with SI BC was already suspected by David 

but as his argument was formulated in the framework of (renormalized) 
continuum PT his claims were open to a counterattack and the issue 
remained unsettled. Our result can be viewed as the lattice analogue of 
David's statement. 

Having identified the SI BC as "sick" , the next question is whether the 
conventional BC are "healthy". Consider the bare lattice PT expansion 
of some short-distance observable on a finite lattice |A| ~ with various 
BC's. Let c^(xi, . . . , L), r > 1 be the coefficient of in an asymptotic 
expansion with BC of type a, where the sites xi, . . . are fixed and a 
distance 0(L) away from the boundary. 

Q2: Which BC's a can be anticipated to give finite and coinciding answers 
for the L ^ oo limits of their perturbative coefficients c"(xi, . . . , x^; L)? 

Q3: Do the infinite volume short-distance observables admit an asymptotic 
expansion in If "yes", can their coefficients c'^{xi, . . . , r > 1 
be obtained from PT via limj;,^oo c"(xi, . . . , Xn',L) = c'^{xi, . . . , 
where a is one of the BC meeting the conditions in Q2? 

In section 3 we shall consider these deeper questions raised in ^] and argue 
in favor of the following picture: (i) All BC's involving only ferromagnetic 
couplings and leaving the interior spins unconstrained meet the condition in 
Q2 - provided free and Dirichlet BC meet it. (ii) Assuming the latter also 
the answer to both parts of Q3 is affirmative. We do not have a rigorous 
theorem but our considerations may constitute a strategy for a future proof 
thereof. Section 4 is devoted to some further discussion. 



We have verified, however, that if we compute quantities well away from the bound- 
aries and the center then the same results are obtained as for conventional BC. 
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2 IR divergence of PT with superinstanton BC 



We consider a two-dimensional square lattice A = {{xi,X2); —L/2 + 1 < 
Xfc<L/2 — l},La positive even integer. The set of points surrounding the 
square, I = L/2, \x2\ < L/2}[j{\xi\ < L/2, |x2| = L/2} will be referred 
to as the boundary dA of A. The 0(n) spin field S^^ defined on A U dA is 
an n-component field of unit length = 1. We restrict attention to the 
standard lattice action 

^(5) = /3^(l-S,-S,+^). (2.1) 

For fixed L the perturbative expansion of the two-point function is 

C"(x, y; P, L) ^ (S, • S,) = 1 - ^ I3~'c'^{x, y- L) , (2.2) 

r>l 

where the superscript a indicates the dependence on the boundary condi- 
tions under consideration. In this section we will only consider two sorts of 
BC's, the Dirichlet boundary condition = 5an for x G dA and the so- 
called superinstanton (SI) boundary conditions [Q]. The latter are Dirichlet 
BC's with the additional constraint that the field at some point zq is also 
held fixed parallel to the fields at the boundary. 

The purpose of this section is to show, by explicit computation, that PT 
with SI BC's is infrared (IR) divergent at third order, while PT with Dirich- 
let or periodic BC is IR finite at the same order. For completeness, and to 
introduce the notations, we first reconsider the first and second orders in 
some detail. Bare perturbation theory can be performed by parametriz- 
ing the fields byQ 5" = /3~2 7r", a = 1,...,A^ where N = n — I, and 
5" = (1 — /3~^7f^.)2. The measure Hx^KdSx^i^l^ — 1) is then replaced by 
exp(-^measure)nj,gA(i7ri, with ^measure = \Y.x ln(l - /3"^vf^). The Feynman 
diagrams contributing to the two-point function are the same for the two 
BC's, only the expression for the free propagator differs. We denote the free 
Dirichlet propagator by D{x,y) and the free SI propagator by G{x,y). One 
has 

0(..i,)=D(x.ri- °'-„7'°'';''" . (2.3) 

^{ZQ, Zq) 

^This is not indispensable, a parametrization independent definition could be given via 
the Schwinger-Dyson equations. 
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In the following we will take the point zq to be the origin zq = (0, 0). Also 
we denote by zi the point (1,0). One easily sees that each c'^{x,y;L) is a 
polynomial in iV = n — 1 of degree r or less, and we write 

< = E<;.-^'- (2-4) 



2.1 First order 

In first order we have, for the Dirichlet BC, 

ci''{x, y) = \n [d{x, x) + D{y, y) - 2D{x, y)] . (2.5) 

As mentioned above, for the SI BC one simply replaces D by G. For the 
special case of the energy expectation value we have 

1 1 

where 

Ao = D{zo, zo) - D{zi,zi). 

cP"' approaches its asymptotic value with power corrections since Aq 
0(1/L^) for L — > cxD. On the other hand 



c?''^(zo,^i) = iv(i-iAo), (2.6) 



cf(zo,^i) = A^^G(2i,^i) (2.7) 

approaches the same asymptotic value extremely slowly with corrections 
1/ In L since 

G(zi,zi) = ^-Ao-^Z-S (2.8) 

with 

Z = D{zo,zo) ~ (27r)-^ InL , for L ^ oo . (2.9) 



2.2 Second order 

At the next order there are three contributions 

c?(^,y) = Ecf M(x,y) (2.10) 

r=l 
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X- 



Figure 1: One-loop diagrams contributing to the spin-spin correlation func- 
tion 



corresponding to the diagrams depicted in Fig. 1 (diagram 2 comes from the 
measure part of the action). One has 

1 



-N 



+ 



D{x,x) - D{y,y) 
iiV [d{x, xf + D{y, yf - 2D{x, yf 



cf^^^){x,y)= -Nj2D{x,y; 



(2.11a) 
(2.11b) 



„Dir (3) 



J2 {d{x, y; i)[D{i, i)D{x, y; i) - D{iJ)D{. 



x,y;j) 



{hi) 



+ 



^ j)} , 



(2.11c) 



where the sums in the last equation are restricted to the case when i and j 
are nearest neighbors and we have introduced the notation 

D{x, y; i) = D{x, i) - D{y, i). 

Let us just inspect the coefficient of N'^. Here diagram 2 does not contribute 
and for diagram 1 one has simply 

cf,2 ^'Hzo,zi) - «(zo,zi) = l[G{zi,zif - Al] (2.12) 
which by eq. (|2.^) tends to 1/32 as L — > oo. For diagram 3 we find 



"SI (3), 



05:2 ^">{zo,zi)-c^.^ (^)(zo,^i) = --Ao + hZ-' + t2Z--' + t3Z-' , (2.13) 

lb 



1 
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with 



where 



^3 = Dizo,if E{i), (2.14) 



\2 



E{i) = Yl [D{i,zo)^-D{j,zo) 
F{i) = [D{i,i)-D{j,j)] . (2.15) 

Here the sums are taken over sites j which are nearest neighbors to i. The 
decomposition in the form of eq. (2.13) is made in order to be able to extract 
the asymptotic behavior reliably. Our ansatz (which we have not proven 
analytically) is that the functions ti have an expansion of the form 

R ■ 

ti{L) ~ Y ^^[(2^)"' InL]^' + 0(1/L), (2.16) 

r=0 

with Ri finite. To obtain the leading behavior we have computed ti{L) over 
a large range of L (up to L = 220, using extended precision arithmetics) 
and taken logarithmic derivatives with respect to L. Our findings are that 

i?i = i?2 = , (2.17) 

and hence the contributions in eq. ( |2.13| ) involving ti,t2 vanish in the limit 
L — s- oo. However, 

One can in fact understand the result for ^3 analytically in the following way. 
Far from the origin and the boundary, i.e., for 1 <^ \i\ <^ L the Dirichlet 
propagator is well approximated by its continuum behavior 

D(zo,i) ~ (27r)-Mn-^, (2.19) 

H\ 
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and so 

E{i) ~ -a^,){D{zo,if) ~ -2(27r)-2|i|-2, (2.20) 
where Of^--^ is the laplacian in the variable i. Thus for we have 



— {2tt)-^ /'"''"^ d^xl^r^ln^ \x\ = - — {2TT)-^ln^L. (2.21) 
32 J\x\>i/L 96 



This is the same result as was guessed in [^. These authors were quite 
fortunate to get the correct result because initially they just evaluated the 
full function C2.2 numerically and then made a rather naive extrapolation. 
We stress that the analysis has to be done extremely carefully on the lines 
outlined above to properly treat the appearance of polynomials of logarithms 
in the denominators 

The final result of this subsection is that 

hm [cll^izo, z,) - c°|(zo, zi)] = ^ (2.22) 

(we also know Huil^oo c^-2{^o, zi) = 0). The Mermin- Wagner theorem im- 
plies that the perturbative expansion of the two-point function is indepen- 
dent of the boundary conditions in the infinite volume limit. Accepting this 
fact we must conclude that the interchange of the limits /? ^ 00 and L ^ 00 
is not permissible for at least one of the boundary conditions involved. We 
now go on to show that this is indeed the case for the SI BC because one 
encounters an infrared divergence. 



2.3 Leading contribution to the third order 

For the purpose of showing an IR divergence at some order r, for generic n, 
it is sufficient to show that one of the coefficients of this polynomial in N 
is IR divergent. We claim that this is the case for the coefficient of A^'^ in 
c^{x,y;L) with the SI BC. 

■^Indeed, 02^2 ^'^^ is not monotonic in L and has a minimum at a large value of L ~ 120. 
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Figure 2: Two-loop diagrams contributing to the spin-spin correlation func- 
tion with terms cx 



Out of 12 diagrams only the 5 shown in Fig. 2 have a nonvanishing 
0{N'^) coefficient. Explicitly 

c|3(x,y;L)=^C3%M(x,y;L), (2.23) 

r=l 



with 



4^3 = Y^iGix,x) - G{y,y)f [G{x,x) + G{y,y)] , (2.24a) 
ci;3 ^'^ = -l[G{x,x)-G{y,y)]J2[G{i,i)-G{j,j)] 

X [G{x,if - G{y,if - G{x,jf + G{y,jf\ , (2.24b) 
4-3^'^= I E {[G{i,i)-G{jJ)][G{k,k)-G{l,l)] 

{i,j){k,l) ^ 

X G{i, k) G{x, y; i)G{x, y;k)+ 3 perms | , (2.24c) 

4h^'^= I E {[Gix,y;zf-Gix,y;j)'] 

{i,j){k,l} ^ 
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X [G{k, k) - G{1, 1)] G{i, kf + 3 perms | , (2.24d) 



+ (i ^ j) . (2.24e) 



Here "+3 perms" stands for three similar terms with summation indices 
permuted according to ^ j) + {k <^ I) + {i j, <-> /)". 

Again we restrict attention to the energy expectation value at the center 
of the lattice. First we note that for the analogous Dirichlet expression one 
can show limi_^oo 63^3 (zQ) -^i) = 0. For the SI functions we again insert the 
formula, eq. (|2.3| ), for G and separate the coefficients of Z~^. 



Using the results in the last subsection together with the knowledge that 
lim/,^00 J2i^zo ^{^0, zi] if'F{i) = it is easy to show that the contributions 
from diagrams 1 and 2 in Fig. 2 have finite limits as L ^ cxo. It turns out that 
all the remaining diagrams 3-5 in Fig. 2 diverge logarithmically. Consider 
first diagrams 3 and 4. Dropping all terms which have finite L — > 00 limits 
one obtains the decomposition 

c% (20,^1) ~ S4^-^ - 32tiZ-5 , 

where is a one-loop sum 

*4 = ^ E^(^o,i)'i?« , (2.26) 

with 

H{i) = [D{zo,zi-if -b{zo,zr,jf] . (2.27) 

i = (i> 

The Sr are more complicated two-loop sums 

*2 = \ Y.Y.H{t)E{k)D{i,k)D{zo,i)D{zo,k) , 

«4 = E E^W^(^)^(^'^)^(^o,i)i^(^o,A:). (2.28) 

i^ZQ k^ZQ 



12 



We have numerically computed all these sums over a large range of L (up 
to L = 210) and their asymptotic behavior as L — > oo was determined 
using extrapolation techniques described in the previous subsection. In all 
cases we found results consistent with those obtained by substituting the 
behaviors of the functions at large separations: 



S3 ~ 



1 

32' 
1 



-Z 



192 



S2 ~ 
■S4 ~ 



-z' 



1 

96' 
1 

240' 



-Z' 



(2.29) 



We do not give all the derivations but illustrate the procedure again with 
the interesting case of S4. Apart from the approximations ( 2.19| ), ( p.20 ) we 
introduce D(i,k) ~ (l/27r) ln(L/|i — k\) so that 



S4 



32 



1,, In ^ In ^ 

i'^k^ \i\ \k\ \i — k\ 



(2.30) 



Now by symmetry we can consider |A;| < \i\ and averaging first over the 
angle between i and k we obtain 



S4 



.0 fr^ h\ Ji r r 



In 7— In ■ 



L 



16 
1 

32 
1 



if^zo 



l(2vr)-6y ^In^^ 



i^zo 



32^ ' L3 



rl_l 

L3 5 



In^L 



In^L 



240 



\i\ 



(27r)-5ln^L . 



(2.31) 



The corresponding numerical result is shown in Fig. 3. The fifth logarithmic 
derivative of S4 multiplied by (27r)^/5! is shown as a function of 1/L. The 
definition of the lattice logarithmic derivative used is 



L 



/(L + 1)-/(L-1) 



The data are calculated for even values of L in the ranges 8 — 60, 100 — 120, 
200—210 using extended precision arithmetic. The dotted line is the analytic 
result 1/240. 
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0.00428 
0.00426 
0.00424 
0.00422 
0.00420 
0.00418 



0.00416 - ^ « 



0.00414 - 
0.00412 - 



0.00410 I • ' • ' • ' • 1 

0.00 0.02 0.04 0.06 0.08 

1/L 



Figure 3: The 5-th logarithmic derivative of S4 multipHed by (27r)^/5! plotted 
versus 1/L. 

The coefficient C3.3 ^^\zo, zi) involves new sums but can be treated sim- 
ilarly. Summarizing we find 



1 

1440" 

4^'"(20,3i) ~ ^Z. (2.32) 

Thus the divergent contributions do not cancel and one ends up with an 
overall logarithmic divergence cf^^{zo,zi) ~ (7/2880)Z as L — >■ 00. 

This shows that one cannot interchange the limits (3 ^ 00 and L ^ 00 
for the SI BC because the perturbative coefficients diverge at two-loop or- 
der. In particular, the SI BC's do not belong to the class of "healthy" 
ones alluded to in Q2 and Q3 cannot even be addressed. Initially however, 
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the SI example was probably thought to discredit PT on internal grounds, 
stemming from the apparent ambiguities in the PT answers for the infinite 
lattice quantities!^. As such, we think, it is no longer valid. Of course one 
could still argue that the divergent SI answer might even hit a point in sig- 
naling that the quantities in question do not have an asymptotic expansion 
in /3~^. However this possibility is not an argument by itself and has to 
be supported by external means, see, e.g., [||. In the next section we shall 
leave this example and try to develop criteria of what healthy EC's should 
look like and under what conditions an affirmative answer to Q3 should be 
expected. 



3 Addressing the questions Q2 and Q3 

We start by comparing the correlator between two spins for different bound- 
ary conditions. The discussion will be based on the validity of the following 
conjecture: For a spin system with ferromagnetic couplings the correlation 
function C{x, y) = (S^ • Sy) increases if any of the ferromagnetic couplings is 
increased. This conjecture is physically rather intuitive, nevertheless it has 
not been proven yet Q. Consider a spin system with a given ferromagnetic 
interaction on: 

1. a finite lattice A with free boundary conditions, 

2. the same lattice with Dirichlet BC, and 

3. an infinite lattice. 

Validity of the above conjecture implies the following inequalities between 
the corresponding correlators: 

Cf-^(x,y;/3,A) < < C°-(x, y; /?, A) , (3.1) 

for any x, y G A. More generally, the conjecture also implies that on any 
lattice A* containing A U (9A the correlator C*{x,y; P, A*) satisfies 

C''^''ix,y;P,A) < C7*(x, y; /?, A*) < CD''^(x, y; /?, A) . (3.2) 
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From here it follows that C^^°°{x^y\ (3^ A*) increases while C^''^(x, y; /?, A*) 
decreases with increasing size |A*|. Further the systems with periodic BC 



on A* will satisfy eq. (3.2). Note that the SI BC with € A (where zq is 
the site where the extra spin is fixed) does not belong to the class considered 
in eq. (|3.2D, the considerations apply, however, when zq is outside A. 



In the following we shall argue that if the correlation functions of the 
systems 1 and 2 have well-behaved perturbative expansions - in a sense to 
be specified below - then 



(i) the infinite system 3 has an asymptotic expansion with coefficients 
given by the formal L ^ oo limit of the perturbative coefficients of 
systems 1 and 2, and 

(ii) all systems considered under eq. (^]^) have similarly well-behaved per- 
turbative expansions. 

For this purpose we have to make some natural and verifiable assump- 
tions on the perturbative series with free and Dirichlet BC's. To simplify the 
notation we shall return to the square lattice of section 2 and write C"(/?, L) 
for C°'{x, y; /3, A). The distance of the fixed sites x, y from the boundary is 
taken to be 0(L), so that C"(/3,-L) qualifies as a short-distance observable, 
as defined in the Introduction. 

Let us elucidate on what we mean by a "well-behaved" perturbative 
expansion. For fixed L and given boundary conditions a = free or Dirichlet, 
one can calculate C" in PT as a power series in 1//? to some arbitrary order 
k: 

k 

C^{I5, L) = l-J2 C?{L)P-' + RUP, L) . (3.3) 

r=l 

There is perhaps no doubt that for fixed L these are asymptotic series, i.e. 

i?^(/?,L) = o(/3-'). (3.4) 

Consider now the formal L — > oo limit of the coefficients c°(L). On the 
grounds of (the lattice counterpart of) David's analysis |^] one expects that 
for the case of the free, Dirichlet (and periodic) BC's there are no infrared 
singularities and that the convergence to the L — > oo limit is of 0(1/L) (up 
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to logarithms, i.e. 0(1/L • In'^L)). We shall assume that the perturbative 
coefficients obtained actually coincide in the infinite volume limit: 

c"(L) = Cr + O (^j^ ) = free, Dirichlet, (3.5) 

provided the distance of the fixed sites x, y from the boundary is 0{L). One 
has good reason to trust ( |3.5[ ). The free propagator for these and other 
conventional BC has a size dependence of 0(1/L^) in contrast to the SI case 
with an 0(1/ In L) dependence. Thus conventional boundary conditions 
should not produce such peculiar behavior as observed in the one-loop order 
for the SI BC's. Assumption (|3.5| ) can in principle be verified order by 
order. It could also be relaxed |Q (see note added), but we expect eq. ( p.5| ) 
to describe the form of the actual finite size corrections. For r = 1, 2 we have 
verified this numerically. Note that r = 2 is already the source of troubles 
in the SI case. 

Clearly, the problem of establishing the claims (i) and (ii) amounts to 
controlling the remainders in an ansatz for an asymptotic expansion of 
C°°{/3) or for the correlators in eq. (^I^). Let us first consider claim (i) 
and define 

i2r(/3)=C~(/3)-l + ^c,/3-^ (3.6) 

This definition is chosen such that claim (i) is equivalent to the statement 
^k^iP) — o(/3~'^), k > 1. That is C°°(/3) has an asymptotic expansion in 
powers of and its coefficients coincide with c^. For later use let us also 
introduce ^ 

AUP,L) = j2[c?{L)-Cr]p-' (3.7) 

r=l 

and note that one can now rewrite ( |3.1[ ) in the form 

At^^iP, L) + 4-^(/?, L) < R^iP) < Af^iP, L) + R^^iP, L) . (3.8) 



The naive strategy to establish (i) would be to first try to show that 
\R^{P^L)\, a = free, Dir are bounded by some L-independent function of 



o{P ) and then conclude from eq. (3^) that R'^{P) = o{P ). Our main 



technical observation is that one can reach the same conclusion using a 
bound on |i?^(/?, L)[ which is much less stringent, and which needs to be 
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established only in some part of the plane. In detail, we assume that 
the bound 

\R-iP,L)\<B,^^^^^^^, (3.9) 

with some finite and p{k) > has been established in the following region 
of the L-P plane 

L<Li{(3), f3>po where (lnLi(/3))P(^V/3 ^ as /3 ^ oo (3.10) 

and Pq is some (L-independent) constant. Convenient choices for the func- 
tion Li{P) will be given below. Similarly define another region of the L-/3 
plane by 

LoiP) <L, P> Po where P^-^/Lq{P) . (3.11) 

Choosing now any unbounded path L{P) in the intersection of both regions 
( |3.1C1| ) and ( p.llj ) the desired conclusion can be reached immediately: One 
has ^fc(/3, L{P)) = o{P~^) by eq. (|]5|) and the condition on Lq{P). Similarly 
from eq. ( ^^ ) and the condition on Li{P) it follows that R'^{P, L{P)) = 
o[P~^). Then the inequalities (^) imply 

Rf{P)=o[p'^) , (3.12) 

which is basically the required result. Strictly speaking, one should for 
a fixed k refer to the numbers Cr, r < k as the (unique) candidates for 
the coefficients in an asymptotic expansion of C°°{P) in /3~^. Only after 
eq. (|3.9|) and hence eq. (|3.12|) has been shown for all k, the existence of the 
asymptotic expansion follows and the numbers can properly be referred to 
as the coefficents of this expansion, so that claim (i) follows. It is also easy 
to see that starting with eq. ( |3.2p the same strategy applies and eq. ( |3.9| ) 
also implies claim (ii). 

Let us add a number of comments. First one should check that choices 
for Lq{P) and Li{P) exist such that the intersection of the regions ( |3.10|) 
and (|3.11| ) is nonempty for each sufficiently large p. A simple example is 
Lo{P) = Li{P) = P^, which is moreover the choice with about the minimal 
growth that allows one to draw conclusions about the first k coefficients of 
the expansion. Of course any choice for Lq{P) and Li(P) growing faster 
than P^ (and such that the intersection of the regions ( ^.10 ) and (3.11) is 



nonempty) is also sufficient. For example a region containing the path 

L(/3) = exp (in^ (3.13) 
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would be convenient, because it allows one to cover all k simultaneously in 
the above argument. However one should keep in mind that the faster the 
upper boundary Li{(3) grows, the more difficult it may be to establish the 
bound ( |3.8| ) within the region ( 3.10| ). 



A second comment concerns the relation to David's theorem. According 
to (the lattice counterpart) of this result one expects that the infrared diver- 
gent terms cancel in the perturbative coefficients. However, to directly prove 
such cancellations in the remainder R'^((3,L) (e.g., p{k) = in eq. ( p. 91 )) is 
probably a more difficult task. The point of using the correlation inequali- 
ties is twofold: First, the claim (i) can be established without assuming the 
IR finiteness of the remainders for fixed (3, as described before. And second, 
once claim (i) is established, the IR finiteness of the remainders actually 
follows from the correlation inequality (|3.2| ). To see the latter let L{(5) be 
an unbounded path in the intersection of the regions (|3lO| ), (|3lll ) and let 



(/?,L') be any point in the region {(/3,L") |/3 > Po, L" > L{P)}. Then by 
(H) 

< A^''{f3,L') + Rf'{(3,L') < ^D-(/3,L(/3)) +i?D-(/?,L(/3)) .(3.14) 

Since, given eq. (|3.9| ), we know the behavior o(P~^) for the terms with L{P) 
and by eq. also for A%{I3,L'), it follows that R'^{(3,L') = o{l3-^) is 

valid "uniformly" in the region L' > L[f3). (In other words, the statement is 
true for arbitrary path L'{0) > L{(3) with a bound on R'^{(3 , L' [fi)) which is 
independent on how fast L'{0) grows. Note that this is not uniformity in the 
usual sense where one would require a region V > Lq with Lq independent 
of /?.) In a sense this supplements (the expected lattice counterpart of) 
David's argument on the IR finiteness of the FT coefficients. 



It remains to establish eq. ( |3.9| ). Since we have not managed to prove 
this, the proof of (i) and (ii) remains incomplete. However, we find the 
strategy outlined promising and for the rest of this section we shall give a 
plausibility argument that eq. (^^) indeed holds. To explain the argument 
let us follow the expectation of Fatrascioiu and Seller that FT breaks down 
for large lattice sizes L since the spins at large relative distances become 
strongly decorrelated. Consequently, one can question the validity of the 
bound (|3.9| ) in the whole plane, in particular for /3 fixed, L — > oo. 
Recall, however, that in the strategy we used to establish claim (i) exclusively 
lattice sizes of type ( p. 13 ) enter which grow with /? only in lattice units - the 



19 



physical size of these lattices L{(3)/S,{f3) decreases as /? ^ oo. (Here ^(/3) 
is the correlation length growing exponentially in the standard scenario.) 
Consequently, in our stratgey it is sufficient to establish the bound ( |3.9| ) in 
the region (|3.10| ) of the L-j3 plane. As it is shown below, for such mildly 
growing volumes the spins become increasingly ordered with growing (3 - 
hence the above objection does not apply. Roughly speaking, the argument 
estimates the probability of having fluctuations vf^ > 6"^ at least at one of 
the sites and finds that such fluctuations are exponentially suppressed by a 
factor ~ exp(— 7r5^/5/ In L). 

We begin by considering the effect of constraining the fluctuations in an 
auxiliary system. Let us calculate the probability that in a massless free 
theory the fluctuation exceeds some given threshold 5. We shall use again 
the A U dA square lattice of section 2 and the action 

X,fJ. 

In the case of the Dirichlet BC the field at the boundary is fixed to zero, 
while for the free BC we use the remaining global translational symmetry 
to fix, say 0o = 0(x = 0) = 0. (For this BC the site x = is assumed to 
be left out in the sums and products appearing in the expressions below.) 
Consider now the constrained model described by the partition function 

Z{(3,6)= f X{d^^e-^''^^'f'\ (3.16) 

Here the field values are restricted to \(f)x\ ^ ^- Define the correction term 

R{f5,5) as 

Z{(3,5) = Z{(3)[l-R{(3,5)] , (3.17) 

where Z{f3) = Z{f3, oo) is the corresponding partition function for the un- 
restricted case. The term -R(/3, 5) can be interpreted as the probability to 
have \4)x\ > 5, for at least one of the sites x. It satisfies the inequality: 

1 p POO 



Z(0) . 27ri 7~oo \ a + iO a — iO 
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\- J- r da( """^ - ex (-Rl^ 



XX 



Here D^x = {4'x)- '^^^ function F is given by 

5 



= - / due-" 



+ J2 ' 

and has the properties F{0) = 1, < 1, and F{oo) = 0. Finally, one has 

< RiP,d) < |A|exp ( ) , (3.18) 

where -Dmax = max^; -D^x- In two dimensions Dmax ~ (^tt)-^ InL hence 

< R{p, 6) < exp (^-^^^ ■ (3.19) 

In fact, 5^/2D XX is the minimal value of the action under the condition that 
(px = To show this consider the configuration 

<^y = ^'5 for any y. (3.20) 

-^xx 

Obviously, (px = ^ and it satisfies the lattice equations of motion, A(f)y = 
for any y ^ x. The corresponding action value is given by 

-W^y{A^)y = -^, (3.21) 

y ^^XX 

as stated above. It is also easy to show that the n-point functions in the 
constrained model will differ from those in the unconstrained one also by 
exponentially small terms of the form ( 3.1S| ). 

We now return to the 0(n) a model. Again, for the case of the free EC's 
we use the global 0(n) symmetry to fix Sq = (0, 1). Using the parametriza- 
tion 'S>x = (vr^., ±-y/l — vf^) it is easy to verify the following inequality :Q 

1-S, -Sj, > ^(tt, -7fy)2. (3.22) 



■'Note that for S" > the parametrization S = (</), + 0^ gives an upper bound, 

1 — Sa; • Sy < \{<i>x — which also might be useful. 
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As a consequence of eq. ( 3.22| ), the standard action of the a model is bounded 



from below by the corresponding free field action, ^(S) > .4o(7r) and 

^(S)>-^ for \^,\>5. (3.23) 

From here it follows that the probability to have \tTx\ > S at least at one 
site X is suppressed by the exponential factor given in eqs. ( p.l8| ) and ( 3.19D . 



Once again, the constraint influences the correlation functions only by an 
exponentially small amount as /3 — > oo. As anticipated, increasing the lattice 
size not too fast, for example, as in eq. ( ^.13 ), the correlation functions of the 



constrained and unconstrained systems will differ by correction terms which 
decrease faster than any inverse power of /3. Further the bound (|3.19| ) allows 
one to choose a constraint S which decreases with increasing (3 not too fast, 
say as S{l3f ~ ln3/31nL//3. In this case the bound still vanishes faster than 
any inverse power of /5, while the quantity /Svr^, the leading perturbative 
term in PA{S), goes to zero as /3 ^ oo. 

By comparing the constrained and unconstrained systems one concludes 
that the fluctuations in the unconstrained system are essentially bounded 
as 7?^ < const In L//3. In other words a system of mildly growing size 
(e.g., as in eq. ( |3.13| )) becomes increasingly ordered as /3 — > co. As shown 
before, within our strategy it is sufficient to prove the existence of the bound 
(3.9) for such mildly growing L{I3). By the above reasoning this should be 



a simpler task, but still remains to be done. 



4 Discussion 



Let us first discuss the physical picture behind the arguments of the previous 
section in more detail. In [Q] the so-called superinstanton configurations are 
introduced. In analogy with the free field case one considers the configura- 
tion with minimal action under the condition that the spin at the middle (at 
X = zq) is rotated by an arbitrary angle relative to the spins at the bound- 
ary (fixed by the Dirichlet BC). The action value of these configurations 
is 0(1/ In L) in two dimensions, so that they can be viewed as saturating 
the bound (|3.23| ). They play a crucial role in disordering the system in 
the infinite volume. As the authors note, the fact that the energy of su- 
perinstantons goes to zero as L ^ oo implies that, in an infinite volume, 
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they are present at arbitrarily large (3 and disorder the spins forbidding a 
spontaneous magnetization in two dimensions. With the assumptions made 
we have shown, however, that for establishing the correctness of the pertur- 
bative expansion it is enough to consider a mildly growing size with L{P) 
given, e.g., by eq. ( |3.13| ). Under such circumstances the superinstantons 
with |7fa;| > 5 are exponentially suppressed and the error made by restrict- 
ing the integration region to \tTx\ < 6 is exponentially small. (After the 
standard PT steps - rescaling tTx and expanding the Boltzmann factor with 
only the quadratic part kept in the exponent - the dependence on 6 is again 
exponentially suppressed, as expected.) For the correlation function (Sx-Sy) 
only superinstantons with size smaller than |x — y| contribute significantly, 
larger ones will rotate the two spins simultaneously. This contribution is 
responsible for the leading order bare PT result 

((S.-So)2)~i^ln|x|. 

This also explains why the 0(n) invariant quantities are infrared finite in 
PT, as opposed to the noninvariant ones which diverge as some power of 
InL. However the form of the finite size corrections is incorrectly given by 
PT, eq. (p. 51). Indeed, due to the nonperturbatively generated mass the 
corrections should be exponentially small for L ^ 

In d = 1 dimensions the status of PT is quite different. In one dimension 
Dma.x oc L and the correction term is of the form exp(— c/35^/L). Conse- 
quently, the 0(1/L) corrections in eq. ( |3.5D are not negligible compared to 
higher order terms in even for such large sizes when the perturbative 
expansion breaks down, i.e., L{P) cx /?. This is in accordance with the ob- 



servation by Hasenfratz [10| that the limits /3 — > oo and L ^ oo in d = 1 
are not interchangeable. As pointed out in [^], starting from order (3~^ the 
coefficients are infrared divergent. 

The case of the SI BC in d = 2 is very similar to the d = 1 case with the 
usual BC. Indeed, already the tree level result has a finite size correction 
0(1/ InL) which is comparable to higher order terms even for exponentially 
large L(P) where PT breaks down. The correlation inequalities for this case 
read 

C~(x, y; P) < C'^'^ix, y; P, L) < C^'{x, y; P, L). (4.1) 
The first few terms of the perturbative expansion for the nearest neighbors 
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are given by 

CSi(zo,zi;/3,L) = 1-^ 

N 

while the corresponding piece of C°° 

N 



16 In L 



+ 



In^L 



48 + 96 +° 



InL 



+ ... 



IS 



1 



32^ 



+ 



(4.2) 



(4.3) 



Clearly, eq. ( ^.21 ) is inconsistent with the correlation inequality (4.1) for 
> 1 (i.e., n = + 1 > 2) when the formal L — > cx) limit is taken. In the 
regime lnL(/3) < ln^(/9) oc (3, however, the 1/lnL correction term in the 
tree level contribution restores the inequality. 

What is the origin of the anomalously large finite size correction for the 
SI BC? Consider a general SI BC where the spin S^q is fixed to a direction 
which is not necessarily parallel to the spins on the boundary dA. Denoting 
the corresponding expectation value by (0)8^ ) the following relation holds 



(O) 



1 



Dir 



(4.4) 



where F{Szo) the free energy of a superinstanton and Z{P) is the partition 
function for Dirichlet BCj^ In the lowest order in 1//3 and for Itt^qI <C 1 the 
free energy is up to a (for our purposes irrelevant) additive constant given 
by (cf. eq. ( pi]) ) 

nS.o) = TT7^7^vf2 . (4.5) 



2D{zo, zo 



The SI solution gives (cf. eq. ( p.20D ) 



1 



Then the integration over the SI solutions (still without the contribution 
from the fluctuations) introduces a nontrivial expectation value for S 
given by 

i_^+or^i. (4.7) 



20 ■ ^Zl 



16/3 InL 



curious point in the terminology is that the SI BC with S^^ — (0, 1) is the one 
which excludes superinstantons centered at zq. Rather the ordinary BC takes care of all 
superinstantons. 
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The analogous correction term in eq. ( [4 .21 ) — which is a result of the fluc- 
tuations for the case when Szq is fixed — compensates the large finite size 
correction in eq. ( [4.7D , as it should, since no such term appears for the 
Dirichlet BC, eq. (|2.6| ). As eq. (|4.5| ) shows, the natural expansion parame- 
ter by integrating over the SI directions is lnL//3, and it is easy to see that 
in generic ©(A'^ + l) models the higher order contributions in eq. 

Q will be 

of the form \/\r? L ■ {XnL/jS)^. As a consequence, a logarithmic divergence 
in order l//3^ is expected also in the perturbative result with the SI BC, in 
agreement with the explicit computation of section 2. Since only the coeffi- 
cient of the leading power in A'^ has been computed one cannot exclude that 
for a specific A^ the divergence is cancelled against one from a subleading 
power in A^. In particular it would be interesting to see what happens in 
the 0(2) model. 

So far we mainly considered bare lattice PT in volumes whose physi- 
cal size goes to zero in the continuum limit, L(/3)/^(/?) — > as /3 — oo. 
However, this was sufficient, together with the correlation inequalities, to 
argue that the standard bare PT provides the correct asymptotic expansion 
for the system in an infinite lattice. The physically more relevant question 
is whether the renormalized PT provides an asymptotic expansion for the 
correlation function at short physical distances, Xphys = eo^, < e <C 1. 
A positive answer to this question is not an automatic consequence of the 
proposed validity of the bare lattice PT. 



Our considerations did not make use of the "integrability" of the 0(n) 
models. Assuming asymptotic freedom one can establish the existence of a 
Yangian algebra of non-local conserved charges and the absence of particle 



production [12, 13 1. The 0(n) symmetry then basically determines the S- 



matrix amplitudes [14|. This bootstrap S-matrix has been tested (at low 
energies) in lattice studies [|^] and used as an input for the thermodynamic 
Bethe ansatz to compute the exact m/A ratio [^]. The results are also 
consistent with the 1/n expansions |17|]. Finally one can use the bootstrap 



S'-matrix as an input for the form factor approach ||l^, 19 1 which provides an 
alternative nonperturbative definition of the theory. For the 0(3) model the 
results in [|^] strongly indicate that the model thus constructed coincides 
with the continuum limit of the lattice theory at least at low energies. At 
intermediate energies the results coincide with renormalized PT and at high 
energies they are consistent with asymptotic freedom. Two non-perturbative 
constants determined exactly in this approach are again consistent with 
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Monte Carlo data |22, 23, 21]. One cannot help feeling that the most nat- 
ural way to reconcile these facts is the conventional wisdom. That is the 
continuum limit of the lattice theory coincides with the model described by 
the bootstrap approach, which is in turn correctly described by an asymp- 
totically free PT at high energies. 
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Note added A concise alternative to our first exposition is also contained 
in a recent note by Sokal In particular he observes that the assumption 
( |3.5[) on the finite size dependence of the coefficients can be relaxed as fol- 
lows: Assume that there exist coefficients {cr}j.^i and powers e(r) > such 
that 

K{L)-Cr\<0 (l-'^'^A , a=free, Dir 



for r = l,...,k. Further replace the condition on Lq{0) in eq. ( ^.1[)| ) by 
P^'"^ Lq{P)~^'^^^ ^ for r = 1, . . . , /c. He assumes that the bound ( ^ ) has 
been established in the whole plane, but this could also be weakened 
along the lines described in section 3. Thus assuming that the bound (|3.9| ) 
has been established for at least one unbounded path contained in the new 



region ( 3.10 ) he shows that C°°{(3) has an asymptotic expansion with the 
coefficients c^. The proof remains essentially the same: One takes an L{(3) 
for which both A'^ and are o{(3~^) under the assumptions stated and 



concludes R'^{(3) = o{(5 ^) from the correlation inequality (3.1). 
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